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ABSTRACT 

We study the cosmic velocity-density relation using the spherical collapse model (SCM) as 
a proxy to non-linear dynamics. Although the dependence of this relation on cosmological 
parameters is known to be weak, we retain the density parameter in SCM equations, in 
order to study the limit fi,n 0. We show that in this regime the considered relation is 
strictly linear, for arbitrary values of the density contrast, on the contrary to some claims in 
the literature. On the other hand, we confirm that for realistic values of Vlra the exact relation 
in the SCM is well approximated by the classic formula of Bernardeau (1992), both for voids 
{5 < 0) and for overdensities up to (5 ^ 2 - 3. Inspired by this fact, we find further analytic 
approximations to the relation for the whole range 5 e [—1, oo). Our formula for voids ac- 
counts for the weak ilm-dependence of their maximal rate of expansion, which for < 1 is 
slightly smaller that 3/2. For positive density contrasts, we find a simple relation 



V ■ V = 3i/o 



0.6 



(1 + ,5)1/6 - (1 + ,5) 



1/2 



that works very well up to the turn-around (i.e. up to (5 < 13.5 for = 0.25 and neglected 
^a)- Having the same second-order expansion as the formula of Bernardeau, it can be re- 
garded as an extension of the latter for higher density contrasts. Moreover, it gives a better fit 
to results of cosmological numerical simulations. 

Key words: methods: analytical - cosmology: theory - dark matter - large-scale structure of 
Universe - instabilities. 
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1 INTRODUCTION 

The gravitational instability is commonly accepted as the process 
of large-scale structure formation in the Universe. According to this 
scenario, structures formed by the growth of small inhomogeneities 
present in the early Universe. Gravitational instability gives rise to 
a coupling between the density and peculiar velocity fields of mat- 
ter. On very large, linear scales, the relation between the peculiar 
velocity v and the density contrast 5 in co-moving coordinates is 



V-v(x) = -J//(n,A)5(x), 



(1) 



where H is the Hubble constant. [For simplicity of notation, we use 
the notation {Q,, A) instead of (Q,m, ^^a)-] The coupling constant, /, 
carries information about the underlying cosmological model and 
is related to the cosmological matter density parameter, fi, and cos- 
mological constant. A, by 



(2) 



( lLahavetalJll99lh . The linear amplitude of peculiar velocities 
is thus sensitive to fi; on the other hand, it is quite insensitive 
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to A. Hence, comparing the observed density and velocity fields 
of galaxies allows one to constrain Q,, or the degenerate combi- 
natio n 13 = 0.°-^ /b m the presence of so called galaxy biasing 
(e.g. Istrauss & Willicklll995l for a review). This is done by ex- 
tracting the density field from all-sky redshift sur veys - such as 
the P oint Source Catalogue Redshift survey (P SCz.lSau nders et 
|200C|) . or the 2MASS Redshift Survey (2MRS, lHuchra et al. 200: 
- and comparing it with the observed velocity field from peculiar 
velocity surveys. The methods for doing this fall into two broad 
categories. One can use Equation l|T), calculating the divergence 
of the observed velocity field and comparing it directly with the 
density field from a redshift survey; this is referred to as a density- 
density comparison. Alternatively, one can use the integral form of 
Equation l|Tll to calculate the predicted velocity field from a red- 
shift survey, and compare the result with the measured peculiar ve- 
locity field; this is called a velocity-velocity comparison. Velocity- 
velocity comparisons are generally regarded as more reliable, since 
they involve manipulation of the denser and more homogeneous 
redshift catalogue data, while density-density comparisons require 
manipulation of the noisier and sparser velocity data. In both cases, 
the density and velocity fields need to be smoothed in order to 
reduce errors and shot noise. Velocity-velocity comparisons re- 
quire a smaller size of smoothing, of a few Mpc. For example, 
IWillicket al] ( Il997h used a smoothing scale of 3 h^^ Mpc. Such 



2 Bilicki & Chodorowski 



scales are called mildly non-linear, the variance of the density field 
smoothed over the scale of a few Mpc is of order unity. 

Mildly non-linear extensions of Equation ^ have been de- 
veloped by a number of workers. These extensions have been 
based ei ther on various analytic al approximations of non-line ar dy- 
namics teegos & Gellei"l989'; Bemardea u 1992, hereafter IB92|: 
Catelanetal. 1995; Chodorowski 1997.: IChodorowski & Lokaj 
19971 : IChodorowski et alj 119981), or numerical (eithe r N-body 



or hydrodynamic) si mulations jMancinelli et al] 19931; Kudlicki 
et al.2000, hereafter iKaCPeROOh. or bo th jNusser et al.1 Il99ll ; 
I Grama milll993l ; iMancinelli & Yahillll995l ; Bernardeau et al.l999, 
hereafter lB99n . Unlike the linear case ([T), the non-linear relation 
between the velocity divergence and the density contrast at a given 
point is non-deterministic (though in the non-linear regime the two 
fields remain highly correlated). Therefore, for a full description of 
the relation, the conditional means (mean V • v given S and vice 
versa) are not sufficient: one has to describe the full bivariate dis- 
tribution function for V • v and 5, or at least the conditional scat- 
ter. These aspects of the v elocity-density relation were studied by 
IChodorowski et al] j 19981) and more extensively bv lB99l . However, 
in practical applications the intrinsic scatter in the velocity-density 
relation is much smaller than the one induced by observational er- 
rors, Mid the conditional m eans are sufficient. 

IB9^ and lKaCPeROd found that very good fits to the mean re- 
lations, obtained for the mildly non-linear fields extracted from nu- 
merical simulations, were given by modifications of the formula of 
IB921 This formula describes a non-linear relation between initially 
Gaussian, random fields ofV-v and 5, under the assumption of a 
vanishing variance of the density field (so the relation has no scat- 
ter). B92 claimed his relation to be the same as the one exhibited 
in the spherical collapse model (hereafter SCM). In practical ap- 
plications (namely with non-zero variance of the density field), he 
predicted his formula to work well in voids, but 'to become very 
inaccurate for 5 larger than 1 or 2'. 

In this paper we study the velocity-density relation in the 
SCM. The reason for such an approach is twofold. First, to derive 
his formula, |B92 used quite sophisticated methods (summing up 
first non-vanishing contributions from the reduced part of all-order 
joint moments ofV ■ v and S). On the other hand, the dynamics of 
th e SCM is very simple and should allow to re-derive the formula 
of IB92| in a straightforward way. More importantly, in the SCM the 
relation can be easily extended to higher values of S, with the hope 
that this modifi cation will fit better the results of numerical exper- 
iments of IB99| andLKaCPeROO. The SCM is in principle insensi- 
tive to the variance of the density field (and the resulting velocity- 
density relation is deterministic), but in practice the variance of the 
smoothed density field dictates how high density contrasts can be 
reached. 

The non-linear relation between 5 and /~^V • v (note the 
scali ng f~^) depends very weakly on cosmological parameters. 
IB92| analysed the Jl-dependence of the scaled velocity-density 
relation in the limit (5^) —> and found it to be very weak. 
iBouchet et al.l 0995) showed that second and third order expan- 
sions for 5 and /~^V ■ v depend extreme ly weakly on SI and A. 
IScoccimarro, Couchman & FriemanI i ll 9990 demonstrated that this 
is the case for all orders. Specifically, they showed that perturbative 
solutions for the density contrast for arbitrary cosmology are, with 
a good accuracy, separable: 5„ — D^{t) en(x), where D{t) is the 
linear growing mode for this cosmology and is the spatial part 
of the 71-th order solution for the Einstein-de Sitter model. Using 
the continuity equation one can then prove, by induction, that the 
velocity divergence depends on Q. and A practically only through 



the f actor f(Q ,, A). Most generally, Nusser & Colberg (1998) (here- 
after lNuCo98l) showed the equations of motion of the cosmic pres- 
sureless fluid to be 'almost independent' of cosmological parame- 
ters. The weak dependence of the scaled velocity-density relation 
on the background cosmological mod el has been also confirm ed by 
N-body numerical simulations ( Mancinelli et aljri993l ; lB9^ . 

However, the il-dependence of the equations of motion of the 
cosmic dust stops to be weak when SI ^ 1 (see eqs. 13-14 of 
,,NuCo 98). This regime of SI is not physically relevant, since the 
currently preferred value of Q, is much higher. Still. lB92l derived his 
formula applying the limit SI ^ 0. Therefore, in the present paper 
we will neglect A (setting A = 0), but will retain the S7-dependence 
of the equations of the spherical collapse and in particular examine 
the limit of small Q.. 

The paper is organised as follows. Section |2| presents gen- 
eral assumptions, terminology and basic formulae of the spherical 
model. In SectionjSjwe focus on the factor /, appearing in Eq. (Tjl 
and commonly approximated by Formula or its simplified ver- 
sion / ~ QP'^ . Section |4| contains an analysis of the regime of 
very small SI and presents the resulting universal velocity-density 
relation. In Sections [5] and [6] basing on analytical considerations, 
we derive approximations for the relation between the velocity di- 
vergence and the density contrast respectively for spherical voids 
and overdensities, for realistic values of Q.. These approximations 
constitute the main results of this paper. SectionjTjgives a compar- 
ison of our fits with results of numerical simulations. We conclude 
in Section jS] 



2 COSMOLOGICAL SPHERICAL MODEL 

Let us consider an open Friedman world model (i.e. with SI < 1) 
without the cosmological constant, A = 0. We introduce the con- 
formal time T] related to the cosmic time t by the equation 

c&t 



dt] = 



Ro a 



(3) 



where _Ro = c/{Ho\/^ ~ ^o) is the curvature radius of the uni- 
verse, c and a are respectively the velocity of light and the scale 
factor; subscripts '0' and (used later) 'i' refer to the present day 
and to some adequately chosen initial moment, respectively. Now, 
the time evolution of the scale factor ca n be expressed in terms of 
the following parametric equations (e.g. IPeeblesll 1 980l) : 

a(j7) = A(cosh7; - 1) , t{ri) = B {sinhrj - -q) (?? > 0) , (4) 

where A and B are constants. Moreover, in this model the confor- 
mal time r] is unambiguously related to the density parameter S7 
by 

2 



(5) 



1 + cosh T] 

If we now consider a top-hat spherical perturbation (a sphere 
of homogeneous density embedded in a Friedman universe), it can 
be analysed as a 'u niverse of its own' (as was noted for the first time 
bv lLemaitrelll93ll) with a scale factor ttp which is the radius of the 
perturbation. Introducing the density contrast of the perturbation 
relative to the background, 5, as 



Pp - Ph 
Pb 



(6) 



we obtain two cases to be taken into account. Using the same ter- 
minology for spherical perturbations as for analogous Friedman 
world models, an open perturbation is such that its initial density 
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contrast 5i is smaller than the critical density contrast 5c (the den- 
sity contrast of an Einstein-de Sitter type of perturbation, i.e. with 
f^(p) _ given by: 



1 



(7) 



It can be checked that the density parameter of thus defined open 
perturbation is fi'^' < 1, as expected. These results are valid under 
the assumption that the initial density of the background is suffi- 
ciently close to the critical density (Q,i ~ 1). The factor 3/5 in 
Eq. ([T} comes from the decomposition of the density field into two 
components, one related to the growing mode and the other to the 
decaying one; we assume here the perturbation to be purely in the 
growing mode. For details see Peebles ( 1980). 

The evolution of such a spherical perturbation is governed by 
equations analogous to {4}: 



ap{(j)) = Ap(cosh<^ — 1), t{(l>) ~ -Bp(sinh<^ - 



> 0).(8) 



The normalization factors are such that {Ap/Af = {Bp/Bf. Of 
course, time t is the same for the background as for the perturba- 
tion, which leads to the relation between and r;: 

sinh(^ — <^ = (1 — r)''''^(sinh»7 — 77) , (9) 

where we have used r = Si/Sc (for a detailed derivation see lPeeblei 
1 1980. ). If (5 > then r > 0, so (f> < ri, and vice versa for negative 
6. 

In order to obtain similar relations for a closed perturbation 
(Si > 5c or fitP) > 1), one should make the following substitu- 
tions: 



Ap 



-Ap 



Bp 



iBp 



(10) 



remembering that in such a case Q < (f> < 2n. 

We can now express the density contrast in terms of the pa- 
rameters r; and 4>, using the relation Pp/ph = (a/ap)'^: 



5 = 



/ sinh (j) — f cosh 77 — 1 
\ sinh ri — T] J \ cosh 0—1 



(11) 



for open perturbations and similarly for closed ones, with the use 
of oil: 



/ cosh 



V 



smn r) — f] I \ 1 — cos r 

(cf. iRegos & Gelledll989l : iFosalba & Gaztanagalll998l) . Note that 
always 5 > —1, but in principle the density contrast has no upper 
bound. However, if initially < 5i < 5c, then 5 cannot exceed a 
maximal value which can be calculated taking — > in i ll lb : 



5lin 



2 (coshr?- 1)^ 
9 (sinh ri — rj)^ 



1 , 



(13) 



The above value becomes the minimal value of the density contrast 
for closed perturbations, i.e. it is a boundary value of possible den- 
sity contrasts between closed and open perturbations for a given 
V- 

The linear theory relates the density contrast of a perturbation 
to its peculiar velocity divergence V • v (Eq. [T]!. In the spherical 
model we obtain V-v = 3{Hp — H),v/heie Hp = dp/ap. For con- 
venience we change units and sign, obtaining what will be called in 
this paper the (dimensionless) velocity divergence, 9: 



(14) 



Some simple algebra is sufficient to find the dependence of 
Hp/H on r? and (j), whic h leads to the following expression 
dRegos & Gelleij|l989l : lB99l) : 



sinh (f> (sinh <j) — (j)) f cosh 77 - 



sinh -q (sinh r) — r)) \ cosh </> — 1 



(15) 



valid for open perturbations on an open background; substitution 
(f) ^ i4> gives the relation for closed perturbations: 



= 3 



sin — sin (f>) f cosh r] 



sinh T) (sinh rj — rj) \ 1 — cos < 



(16) 



Both the density contrast and the velocity divergence, as given by 
dl lb and JI5b . or il2i and JI6t . are parametrically dependent on 
(j) (rj is fixed). Our aim here is to eliminate this parameter (at least 
approximately) and to obtain the 9-5 relation in the spherical model 
in an analytic form. 

As a first step, it is useful to simplify the formula for 9 in- 
cluding 'the easy part' of the dependence on 5. This is done by 
calculating (sinh (j) — (^)/(sinh rj — rf) from (TTJ and inserting the 
resultant expression into l |15t . Then, owing to the hyperbolic iden- 
tity cosh^ X — sinh^ x — 1 and the relation lO for Q,, we finally 
obtain a simplified formula for the velocity divergence: 



= 3 



if7(l + 5)(l-Hcosh( 



(17) 



These considerations were valid for open perturbations. If 
ri'*"' > 1, then we have 



1=F A/|f^(l + <5)(l + cos< 



(18) 



where '-' applies to the case Q < cf> < n and to tt < < 2ti. 
Formula ( I17t [(11811 is simpler than U5b rJ16b1. but the dependence 
on (j) remains; the parameter <j) is related to 5 by Equation dl lb 
[01]- 



3 FACTOR / 

The linear theory (valid for small values of 5) relates the veloc- 
ity divergence as defined above to the density contrast through the 
equation 



/(SI, A) is given by 



9 = f5 

[cf. Eq. (HJl], where the factor / ; 

/ = dln-D 
d In a 

The quantity D{t) is the growing mode of the perturbation. 
The fac t or / h as been a subject of s t udy i n many papers (e.g 



(19) 



(20) 



1991 



I99S 



Feeble j Il976l : iLightman & SchechteJ Il99 0': 'Lahav et alj 
Martel| J^199ll |Bouchet et I1995I ; iFosal ba & Gaztafiaga 
NuCo98). The best-known and most widely used approximation 
(often without reference) is the one given by (Peebles, U976,) : 



(21) 



In this part we will compare this fit with the exact formula for /. 

The spherical model as described here allows us to calculate 
f{Vt, A = 0) as the limit 



/ = lim - . 



(22) 
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Figure 1. Factor / = (d In D) / (d In a) as a function of density pai'ameter 
Q for world models with A = 0: exact relation (solid line) and approxima- 
tion / ~ n"'" (dashed line). 



It can be checked that choosing \5\ <C 1 is equivalent to taking 
|r| ^ 1 [Eq. Moreover, from the relation ^ it follows that in 
this case (j) = rj + e, where je| <C 1. Using the first-order approx- 



imation (1 



N3/2 



1 



jr and expanding hyperbolic functions 



around e = 0, we can linearize Equations (|9}, Jl lb and dlSt . As a 
result we get a linear relation between e and r and further on also 
linear dependencies of 5 and ^ on r. Diving thus obtained velocity 
divergence by the density contrast, we get the following formula 
for / as a function of 77: 



3 r; (2 + cosh r?) — 9 sinh 77 



(23) 



3 (cosh Tj + 1) (sinh r] — rj) — 2 sinh 77 (cosh 77 — 1, 

A similar r elation, but for a 'closed' mod el of the background, can 
be found in iLightman & Schechteil(ll99(]l) . If we now make the sub- 
stitution 77 = arcosh (2/57 — 1) [Eq. ([5}], then after some algebra 
we can express the parameter / as a function of fi: 

3^(0 + 2) \n[2n-^{l + 



n) - 1] - isnVT^ 



.(24) 



12^1^ - 8^/(1 - fi)^ - 60 ln[2n-i (1 + ^/T~~n) - 1] 

This is the exact expression for f( Q.) with A = and Q < 1. I t 
was already derived for example by iFosalba & Gaztanagal il998h . 
Figure [T] presents a comparison of this relation with the Peebles' 
formula OP'^. It can be seen that the power-law approximation is 
sufficiently exact, especially for the currently favoured value of the 
density parameter (fi ~ 0.25). Moreover, owing to the complicated 
form of J24t . the latter is not very useful. However, one should al- 
ways bear in mind that the formula ( 121b is merely an approximation 
and in some applications its usage may lead to errors. A much better 
fit is the one given in a footnote of .NuCo98., : / = t74/7+(i-n)='/20 
Its errors relative to the exact value for the model with A = are 
below 0.3% for f7 > 0.1. 



4 LIMIT OF SMALL 

Let us now examine more thoroughly the limit of f2 ^ 0. We begin 
with checking the asymptotic behaviour of f{0.). This regime, al- 
though not physically interesting, allows to take a closer look on the 
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Figure 2. Factor /(H) in the limit H <C 1: exact relation (sohd 
line), power-law approximation (dashed line) and asymptotic relation 
/(Q) ~ -(3/2)S7(lnn + 3 - In 4) -dotted line. Note logarithmic scales 
on both axes. 



bottom-right end of the diagram presented in Fig.[T] and the results 
obtained will be useful later in the paper. (See also Appendix IaI) 
Starting with the relation J23b and remembering that the limit of 
small Q means 77 ^ 1, we get the following approximation: 



fiv) 



-"(7,-3) (77 



(25) 



If we now observe that for such 77 we also have 77 ~ In 4 — In S7 
and e'' ~ 4fl~^, we obtain an asymptotic formula for /(fi): 

3, 



-J7(lnf7 + 3-ln4) (fi < 1). 



(26) 



Figure|2]clearly shows that for sufficiently small Q, i.e. < 0.01, 
the power-law of Peebles could no longer be used. This plot is also 
a confirmation that in some cases the usage of log-log diagrams is 
well-grounded. 

iB92 studied the cosmic statistical relation between the non- 
linear density contrast and the velocity divergence, evolving from 
Gaussian initial conditions, in the limit of a vanishing variance of 
the density field. He found the result 'to be very close to' 



(27) 



Here, and from now on, the so-called scaled velocity divergence, 
O, is defined as 



(28) 



Note that for \ S\ <^ 1, the non-linear formula l |27b correctly reduces 
to Q^^ -* — 6 , i.e. to the linear-theory relation ( |19b . As already men- 
tioned, |B93 claimed his rel ation to be the same as the one exhibited 
in the SCM. In turn, 'B99' argued that the approximation ( |27b 'is 
strictly valid in the limit Q ^ 0' . Here we check these statements, 
applying the regime 57 > to the equations of the SCM. 

If 57 ^ 1 then 77 ^ 1. Therefore, since for voids (5 < 0) 
we have (p > rj, also 1. For overdensities (5 > 0), the limit 
77 — > 00 applied to Eq. l |13b gives Si^-n —> +00. Thus we can focus 
only on Formula dl lb for S. From Eq. i ll lb we see that in order 
to keep 5 finite (though arbitrarily large), also (j) should tend to 
infinity. In other words, if 77 3> 1, then also 0^1, both for voids 
and overdensities. Hence, still from Equation dllb . we get, up to the 
leading order. 
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<t>i^V- + S) , (29) 
and up to the second order 

02 = r? - ln(l + S) + [4(1 + <5) ln(l + S) ~ 5(4r? - 6)] e'^XSO) 
Equivalently, 

(1+5) cosh 02 = cosh»7+i [4(1 + 5) ln(l + 5) - 5(4r? - 6)].(31) 

Applying this formula in Equation U7b and using the large-77 limit 
of the function /(ry) (Eo.l25b we obtain 



e 



N{S) 



5 + 



N{5) 



lnn + 3-In4' 
where 

N{S) = (5 + l)ln(5 + 1) -5. 

In the limit ^ the second term in Equation 
hence 

e = 5. 



(32) 

(33) 
vanishes, 

(34) 



This is exactly the same relation as for the linear regime (where 
\6\ <^ 1). However, here th e de nsity contrast can have an arbitrary 
value. Thus, the formula oflB92 J27t does not describe the dynam- 
ics of perturbations in the limit ^ 0. The relation l |34t . being 
very simple, is a non-trivial result. When Q, tends to 0, then also the 
(not scaled) velocity divergence 6^0 (peculiar velocities van- 
ish with diminishing Q). However, the quantity 0, as introduced 
by Eq. J28b . converges to a non-zero value for — > 0, due to the 
presence of the factor f{0.), approximated by J26b for small f2. 
The normalisation used here leads to = (5 in the linear theory. 
Why for very small values of fl this relation holds also in the non- 
linear regime? It turns out that this is a general result of dynamics 
in a low-density universe, and does not rely on any symmetry. The 
derivation is presented in Appendix IaI 

For large but finite values of rj. Formula i32\ can be applied. 
Specifically, it can be used for 77 significantly greater than 3 (SI sig- 
nificantly smaller than 0.2), which falls well below the presently 
accepted value of the cosmic density parameter. Therefore, the ap- 
proximation J32b is of no practical importance; we have to continue 
our search for a relevant relation. Just for i llustrative purposes, on 
Figure[3]we plot the exact relation, the lB92l approximation, and the 
approximation i32\ , for an exemplary value of SI = 10"'' (r; ~ 13). 
We see that although for this value of the relation is still non- 
linear, the B92 approximation drastically overestimates the degree 
of non-linearity. 



5 RELATIONS FOR VOIDS 

As voids we will understand any underdense perturbations, i.e. 
those for which 5 < 0. In this section we examine the behaviour 
of the velocity divergence vs. the density contrast for such inhomo- 
geneities. 

When considering overdense perturbations (with 5 > 0), the 
regime of 5 ~ 1 is usually called weakly (or at most mildly) non- 
linear. It ma ^ thus s eem that it should be similarly for the limit 
S> -1 (cf.' Martel|[T99lh . However, if we analyse Equation (TTJ, 
which is valid both for voids and for open overdensities, we can 
see that for finite values of (which correspond to non-zero Q,), the 
condition 5^—1 may only be satisfied for (f> — > -|-oo. Hence, 
the evolution of such a perturbation is highly non-linear when the 
density contrast approaches its minimum value. 




Linear theory 0=8 
Exact relation 0(8) 
Smali-n approximation 
Formula of Bernardeau 



1 2 3 

Density contrast 8 

Figure 3. A comparison of the relation between density contrast S and 
scaled velocity divergence © = /^^ 6 for very small Q (here Q = 10^^): 
linear theory result © = 5 - dotted line, exact relation @{S) - solid line, 
approximation for small f2 - crosses and lB92l approximation )27t - dashed 
line. 



The scaled velocity divergence O, as defined in Eq. ([28), is 
a monotonically increasing function of 5 (for rj, or fl, treated as a 
fixed parameter). Its minimum value is Omin = 0(— 1) (dependent 
on rj), obtained easily by calculating the limit +00 in l ll5t : 



Qn 



(cosh f] -— 1)^ 



sinh ri (sinh rj — rj) 



(35) 



For rj ^ 0, equivalent to S7 — » 1 (the Einstein-de Sitter model of 
the universe), we get the value of Qmin ~ —1.5, which can also be 
calculated otherwise^ The opposite limit of S7 ^ (77 ^ +00) 
leads to Omin ^ —1; this can be equally deduced from l l34t . If we 
adopt the currently accepted value of Q p ~ 0.25 {rjo ~ 2.63), we 
obtain 0(— 1) ~ —1.43. Thus, the lB92l approximation l l27t . which 
gives 0(— 1) — —1.5 independently of SI, has a relative error of 
approx. 5% in this limit for such Qq- 

We would now like to fi nd an (approximate) relation 0~5 for 
the whole range S G [—1,0]. B92 derived his formula expanding 
the relation around S — 0. We adopt a different approach: we ex- 
pand the relation around 5 — —1. (That is, at a first step we assume 
< 5 + 1 ^ 1). Then, for arbitrary rj, the perturbation parameter 
0^1. From Equation i ll U we obtain coshi^i = fl('7)/(l + 
where 



giv) = 



(cosh r; - 1)^ 
(sinh 77 — 77)^ 



and further on 

cosh (j)2 = cosh 01 + 3 — 201 



(36) 



(37) 



^ In the E-dS model we have H = 2/{3t) and /(H) = 1; adopting the 
empty world model (Milne model) for the perturbation, we get Hp = 
and further on e = 9 = 3(1- HpH~^) = -1.5. 
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Using Equation l l37t in Equation l ll7t we get 



62 = e„ 



3(sinhr?-77)(0i -2) 



sinh ri (cosh — 1) 
wliere 

01 =ln[23(r,)]-In(l + 5). 



(1 + 5), 



(38) 



(39) 



Equation ( 1381 ) satisfies explicitly the highly non-linear limit 
G(— 1) = Smin- Also, in the limit 7; 3> 1, this Equation reduces to 
asymptotic Equation J32| >. as expected. 

The range of applicability of formula ([38} is very limited: it 
starts to deviate from the exact relation for 5 about —0.9. We would 
like to introduce such a modification so as to satisfy also the linear- 
theory limit: for \5\ <ti 1,0 = S. Therefore, we adopt the following 
three boundary conditions: 

A. e(-l) = Oniin, 

B. 6(0) = 0, 

c. (de/d(5)u=o = 1. 

Inspired by Equation J38l >. we write 

e = Qmin + ai{ri){l + 5) + a2(r/)(l + S) ln(l + 5) , (40) 

where ai and 02 are arbitrary functions of r/. Imposing the three 
boundary conditions A.-C. on the above formula we find 



e = 5 + (i + e,ni„)iv(5), 

where Smin as a function of Q [cf. ( I42H is 



(41) 



1 - 



2^1"^ - O ln[2fi-i (1 + ^/T~^) 




0.50. 



,0.12-0.06 n 



(42) 



[here /(fi) is given by ( I24H and N{5) has the form of ( 133b . Indeed, 
formula J41t meets all the three boundary conditions: the last two 
are fulfilled since for small 5, N{S) = 5^/2 + . . ., and the first 
one because 7V(— 1) = 1. This simple approximation is robust for 
S close to —1 and around 0; for intermediate values of 5 it slightly 
underestimates the exact value of O (with a maximal relative error 
of 2% for n ~ 0.25). 

Formula ( I4U is probably already sufficiently accurate for 
practical applications. Still, it is of course possible to improve it. In 
order to do this, we expand the exact Q-5 relation around S = —1 
up to third-order in the perturbation parameter (f>. The result is the 
following series: 

03 = 0.0 + ai(l + (5) +02(1 + (5) ln(l + 5) +03(1 + (5)^ 

+ 04(1 + Sf ln(l + 5) + 05(1 + Sf lri'(l + <5), (43) 

where ai are some functions of rj (see Appendix iBt. From the six 
terms above we construct their linear combinations which fulfill the 
constraints A.-C. This, together with the condition of simplicity, 
leads us to postulate 

e = S+[l + QnM^)] N{S) + ai{n) S{l + S) ln(l + S) + 

a2in){l + 5fln''{l + S), (44) 

Fitting this formula to the exact relation gives ai = 0.12 and 
Q2 = —0.09 for Q = 0.25. The fit is very accurate: it has a max- 
imal error smaller than 0.2%. In general, both ai and 02 depend 
weakly on fi: we have ai ~ 0.19 fi"'^" and 02 ^ -0.15 S^"'^^ for 




Linear theory = 5 
Exact 0(5) 
Fit for voids 
Formula of Bernardeau 
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Density contrast 5 



-0.2 



Figure 4. Relation between density contrast and scaled velocity divergence 
for spherical voi ds (u nderdensities): linear theory (dotted line), exact rela- 
tion (solid line), |B92| approximation (dashed line) and fit given by Equa- 
tion I44t (crosses). The density parameter of the background is f2 = 0.25. 



0.1 < n < 0.9. Figure |4] presents a comparison of the exact rela- 
tio n for Q(5), calculated for voids from l ll7b . with the fit ( I44t and 
the B92 approximation ( I27t . As w e can see, our fit lies accurately 
on the approximated curve and the lB92l formula slightly underesti- 
mates exact values of O for 5 close to —1. However, it should be 
admitted that the latter is considerably simpler than ours. 



6 OVERDENSITIES 

An overdensity is any perturbation for which 5 > 0. As already 
mentioned, these can be of two types, depending on the initial 
density contrast: 'open' or 'closed'. For a specific value of Q. (or, 
equally, rf), the boundary value of the density contrast, maximal for 
the first type and minimal for the second, is given by ( 1131 ). For the 
currently accepted value of ilo — 0.25 we have Sn^^ ~ 1.6: such 
overdense but open perturbations {5 < (5iim) fall within the weakly 
non-linear regime. 

In order to find an approximation for 0{S) for overdense 
spherical regions, we will use a similar procedure as we did for 
voids, examining the highly non-linear regime (S 2> 1). Owing to 
the considerations above, it is sufficient to focus on closed pertur- 
bations; the formula for O is then of the form l ll8b . with the '-I-' 
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sign. Highly non-linear infall means that the overdensity collapses 
to a point: the conformal time of the perturbation 4> — > Stt. This 
is in general not physical, as in practice for </> < 27r virialisation 
would occur and prevent further collapse. However, as in the case 
of voids, examination of this regime leads to interesting formulae. 

First of all, we can directly put (p — 2n into dlSl l. getting the 
'1-st order approximation': 



(45) 



We can see that the |B92| formula l l27t . which was not intended to 
work in this regime, indeed will not work: already the slope of the 
curve is incorrect (2/3 instead of 1/2). For realistic values of il, 

{Q,)y/Ti ~ Using this approximate equality and ne- 

glecting the constant term in Equation l |45l l yields the '0-th order 
approximation', Oo = 3 + S. The same relation can be 

deduced from dynamical considerations ( namely, f rom energy con- 
servation in the highly non- linear infall). INuCo98| also found such 
a form of the weak Sl-dependence of the peculiar velocity 

in virialised regions. This is not surprising, since both in our and 
their case, 5 2> 1 and O <C <5. 

Expanding the relation dlSb around (f) — 2n (0 < 2n ~ (p <^ 1) 
to higher order, we obtain the following series: 

e = 3/-'(n) [i + VTi{i + 5)^^^ + ai/6 + S)^^"" 



+ 



+a. 



1/6^ 



+ ■ 



(46) 



where ai are functions of Q only. In order to obtain a fit that would 
both converge to l l45t in the highly non-linear regime of 5 3> 1 
and have proper behaviour in the vicinity of 5 = (conditions B. 
and C. from Section|5]l, we proceed similarly as we did for 5 < 0. 
First, already here we neglect the fourth (and all next) component 
of the series. Then we modify the expansion l |46b by introducing 
two parameters A, B and an integer n: 

e = 3/"' [yi + V^\/rT5-B\/n(l + 5)'''"] . (47) 

Imposing the conditions 0(0) = and (AQ/A5) |i=o = 1> we 
obtain: 



A 



3^ 



f{Q.). (48) 



In particular, for n = 6 [cf. J46l >1 we have 

A{9.,n = &)~2QP''{l~Q.°-^) , 
and 

B{n,n = Q) ~?,-2Q!^-^ 

(remembering that / ~ OP'^). Inserting the above expressions for 
A and B into Equation ( 147 1 with n = & and neglecting the weak 
f2-dependence [since f^^ {yi)\/Q, ~ 57^"'^], we obtain the fol- 
lowing result, which can be treated as a generalisation of the lB92l 
formula: 



(49) 



(50) 



(51) 



An interesting feature of this fit is that it has the same second- order 
Taylor expansion around 5 = as the approximation given bv lB92l : 

e = 5 - i<52 + ■ ■ • . (52) 
5 

This means that in the weakly non-linear regime these two approxi- 
mations give similar results. However, for mildly non-linear values 
of 5 (from 5ii^ up to the turn-arouncfl) our approximation works 

^ The turn-around of a closed perturbation is the moment when it stops 



a 
o 
1) 



' 1 ' 
- £2=0.25 
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Figure 5. A compaiison of proposed approximations for ©(5) with the 
exact relation for spherical overdensities in the mildly non- linear regime 
(up to the turn-around). Solid line shows the exact relation, the|B9^ approx- 
imation j27t is illustrated by dashed line and crosses present the fit given by 
(ST); dotted line is the linear theory relation. The density parameter of the 
background equals to f!o = 0.25; the density contrast of the turn-around is 
then <5ta — 13.5. 



generally better than the formula of lB92l The maximal error of our 
fit for such an interval of density contrasts is about 1.5%. Figure|5] 
shows the discussed approximations for the weakly and mildly non- 
linear regime. 

For higher values of 5, neither the B92 approximation, nor the 
fit l |51| l are adequate. In case of the first one this is mainly due to a 
wrong slope of the curve; in case of the second - due to the neg- 
ligence of the dependence on 57. For that reason in the regime of 
very big 5 we prefer to use the fit l |47t , of a more general form. The 
approximation l |5U suggests that the best choice of n is n = 6; 
however, it turns out that in practice, for highly non-linear density 
contrasts (greater than the value for the turn-around), approxima- 
tion ( 147 1 with n = 4 works slightly better than with n = 6 (with 
the weak S7-dependence included in both cases). In Figure |6] we 
show the behaviour of the function 0(5) in the highly non-linear 
regime. For comparison, we plot the formula of lB92l the simple 
approximation ( I51t and the approximation ( 147 1 with n = 4. 

Our results for the highly non-linear regime are rather of aca- 
demic value, since, as stated earlier, highly non-linear infall is con- 
siderably modified by the effects of virialisation. To account for 
them (and for deviations from spherical symmetry). IShaw & Motal 
(2008) constructed an improved (extended) semi-analytical spheri- 
cal collapse model. For 5 up to about <5ta (which equals to ~ 4.6, as 
the background assumed in the discussed paper is of the Einstein- 
de Sitter type) their model coincides with the standard spherical 
model (studied here), while for larger density contrasts it deviates 
from the latter and under some additional assum ptions matches 
well th e results of N-body simul ations presented bvlHamilton et al.l 
( Il99lh . Indeed, formula (20) of IShaw & Motj ( |2008|) . for T ^ t 
(the limit of the standard model), reduces to our Equation ( I18l l 



expanding, i.e. (ip = 0. In the spherical model as discussed here it oc- 
curs for (/) = vr; the density contrast for the turnaround spans from 
= 97r2/16 - 1 ~ 4.6 for C = 1 to 5ta 30 for H = 0.1. If 
n = 0.25, then 5ta 13.5. 
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Figure 6. An illusti'ation of the behaviour of the function 0(5) for spheri- 
cal overdensities in the highly non-linear regime, i.e. up to the virialisation. 
The solid curve is the exact relation, the dashed line shows the [B92i approx- 
imation and crosses present two approximations described in the text: plus 
marks represent the formula (sT) and multiplication mai'ks show the tit (47) 
with n = 4. 



(their hsc = ^/3). The authors argue that for background uni- 
verses with dark energy their formula is valid only for 5 > 100. 
They claim that for smaller values of 5, their results are not ac- 
curate. We disagree with these statements. As already stated, the 
weak Q. and A dependence of the scaled velocity-densi ty relatio n 
has been shown on the level of the equations of motion 
so independently of the level of non-linearity. Since for small red- 
shift s dark energy beh aves similarly to the cosmological constant 
(e.g. lRiess et al.ll2007h . and since only for such redshifts the weak 
dependence of equations of motion on cosmological parameter s 
starts to play any role (because earlier we had n < 1; lNuCo98h . 
the velocity-density relations for cosmological models with and 
without dark energy must be similar. 



7 COMPARISONS WITH FITS TO NUMERICAL 
SIMULATIONS 

iKaCPeROd studied the mildly non-linear velocity-density relation 
using the Cosmological Pressureless Parabolic Advection (CPPA) 
hydrodynamical code. They found that the mean relation between 
the scaled velocity divergence and the density contrast can be very 
well described by the so-called '7-formula', 



e = 7[(i + 5)'''^-i] 



(53) 



with 7 ~ 1.9. This formula is a modification of the IB92| for- 
mula with 7 instead of 3/2. The offset e > is introduced to 
account for an effect of a finite variance of the density field: the 
value of e is such that the global mean of O is zero, as required. 
(Another effect of a finite variance is to modify the degree of 
non-linearity of the relation.) Without the offset, the above for- 
mula yields 0(— 1) = —7 = —1.9 for 7 = 1.9, in signifi- 
cant difference with the value —1.5, obtained neglecting the weak 
f2-dependence of the exact limit, Eq. l |35b . However, for Gaussian 
smoothing scales of a few Mpc, employed in lKaCPeROOl the offset 
shifts the value of 6(— 1) much closer to —1.5. 



|B99| analysed the velocity-density relation using N-body si- 
mulations performed for various background cosmologies. They 
noticed a weak dependence of the relation on Q, and A. B99 in- 
vented a somewhat more elaborate fit to the extracted mean rela- 
tion, presented in the form of density in terms of velocity diver- 
gence. 



5 = /3(i + e/7r 



1 . 



(54) 



Here, /3, slightly smaller than unity, plays a role of the offset e in 
Equation J53t : it assures that the global mean of 5 is zero, as re- 
quired. In Equation ( I54t 7 is not a constant, but is approximated as 
a following function of S: 



7 = ^+0.3n'''« 



e.2 



(55) 



The above equation quantifies the fact that for larger values of ve- 
locity divergence, the observed relation becomes more non-linear. 
Indeed, 7 grows with growing O (we recall that 7 = 1 corresponds 
to the linear theory). Moreover, for O = —3/2, we have 7 = 3/2 , 
so then 5 = —1, as it was intended. [Note a typo in eq. (20) of|B99|: 
instead of 6 (in our notation, Q), there should be 6.} 

How do these findings, based on fully non-linear simulations, 
relate to our results? In overdensities, our Formula ( 1511 ) fo llows 
closer to the exact relation in the SCM than the formula of IB92| . 
Moreover, our approximation is a formula with increasing e ffect ive 
index 7cft . Its second order expansion is the same as that of B92, so 
for small 5, 7eff = 3/2. For large density contrasts, the second term 
in Equation ( I51t becomes negligible, so asymptotically 7eff = 2 
(for 5 3> 1). T herefore, qualitatively our formula is consistent with 
the fit of lB99l in a sense that 7, as a function of O or 5, is growing. 
It is also consistent with the fit of KaCPeROO, in a sense that the 
average 7 is slightly larger than 3/2. Clearly, our formula is a bette r 
fit to the results of numerical simulations than the formula of lB92l . 

Of course, quantitatively there are discrepancies. First of all, 
it is strictly impossible to satisfy simultaneously the features of 
both fits: 7 is either constant or increasing. This discrepancy be- 
tween the results of the two groups is not necessarily a sign of a 
major flaw in any of their analyses. The two groups used differ- 
ent codes: N-body versus hydro. The first one follows accurately 
non-linear evolution, but provides naturally a mass-, not volume- 
, weighted velocity field, while the latter is needed. CPPA, as any 
hydrodynamical code, provides naturally a volume-weighted veloc- 
ity field, but follows the non-linear evolution after shell crossings 
only approximately. Moreover, the density power spectra used in 
both simulations were different. Also, fit ( l54l of ,B99 was found 
for top-hat smoothed fields, while fit (|53j of lKaCPeROol was elab- 
orated for fields smoothed with a Gaussian filter (more appropri- 
ate for velocity-density comparisons). The effects of smoothing, 
though small, are different for these two filters (see e.g. Table 1 of 
rKaCPeROQ) . Finally, an inverse of the forward relation (density in 
terms of velocity divergence) does not strictly describe the mean 
inverse relation, due to scatter. 

Which results better reflect real non-linear dynamics of 
cosmic random density and velocity fields? Instead of bet- 
ting, it would be probably best to repeat the analysis us- 
ing an output from high-resolution N-body simulations with a 
ACDM power spectrum, employing - instead of a Voronoi tes- 
sellation (Bernardea u & van de W evgaert 1996) - a much sim- 
pler algorithm of extracting volume-wei ghted velocity field of 
IColombi. Chodorowski & Tevssieil ( l2007h . Voronoi tessellations 
are complicated and very CPU-consuming, so they can be ap- 
plied only to a limited number of points, while the method of 
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IColombi et J] ilOOlh can be (and actually has been) applied to all 
simulation points (512'^ in their work). If the actual relation is not 
more non-linear than in the highly non-linear regime of the SCM 
(7 — 2), then we can use Formula J47t . with neglected weak Jl- 
dependence and n treated as a free parameter. Let us write it ex- 
plicitly: 



e 



(2,n) 



i^ + S) 



1/2 



6 



(1 + 5) 



i/ii 



(56) 



where e — n/2 — 3. [For n = 6, it reduces to Formula dSlt.l 
For example, if the best-fit value of 7 is found to be close to 1.9 
and fairly constant, then n = 2.3 would provide an excellent fit. 
Instead, significant 'run' of the index 7 would probably demand 
n > 6. 

If, on the other hand, the results of IB99I are found to be ac- 
curate, then for Q, = 0.25 Equation ( 155b yields 7 = 2 already 
for Q ~ 2.3, and even more for higher G. In this case, in order 
to describe the relation up to the turn-around, one should modify 
also the exponent of the leading term in Formula ( I56t (1/m instead 
of 1/2, with m>2)[f| It is a matter of choice if to fit one 'run- 
ning' exponent (7) or two constant (m and n). In any case, it is 
better to use an additive offset e instead of the factor /3, appearing 
in Equation ( |54t : in applications to velocity-velocity comparisons 
the value of e is not relevant at all. The mildly non-linear veloc- 
ity field is vorticity-free to good accuracy, so the predicted velocity 
field (from the density field) is 



v(r) = 



47r 



(57) 



and the contribution of the offset to velocity averages out to zero. 

An advantage of the 7-formula over Formula l |56| l or its mod- 
ification is that it works also for voids. For underdensities, the for- 
mula of B92 is a very good description of the exact relation in 
the SCM, except for the very tail 5 ~ — 1 (where the weak Q- 
dependence becomes important). Results of numerical simulations 
show very limited need to modify the formula of IB92| for voids - 
the discrepancies appear at larger density contrasts. As stated ear- 
lier, iB9Z predicted t his fa ct. Our formulae for voids give results 
very similar to that of lB92l but describe better the regime 5 ~ — 1. 
This regime is important for predicting expansion velocities of al- 
most completely empty voids (e.g., see TuUy et al. 2008). There- 
fore, using approximation ( I41b for underdensities, we propose the 
following combined formula: 



e = 



Q(m,n) 

5 - 0.5 Q, 



0.12-0.06 n 



S>0, 

N{S) + e, -1<5<0. 



(58) 



Here, Ot"*-"' is the 'm-modification' of Formula i56i , N{S) is 
given by Equation ([33} and e, treated as a free parameter, is the 
same in both cases. Alternatively, as the relation for voids, one 
could use more complicated and extremely accurate Equation J44l (. 
To sum up, we admit that it is disputable if to fit results of nu- 
merical simulations using our formula ( 158b . or 7-formula. What 
remains indisputable is that for overdensities, our standard formula 
(with m = 2 and 71 = 6) is a better starting fit than the standard 
formula of iB92. (with 7 = 3/2). 



^ This modification would create a coefficient of the leading term equal to 
m/2 and modify the offset to e = n/2 — (3m/2). 



8 SUMMARY AND CONCLUSIONS 

The main motivation of this paper was to rederive the formula of 
IB92| in a simple way, using the spherical collapse model (SCM), 
and to extend it to larger density contrasts, where it is no longer 
valid. The undertaken project abounded in surprises: 



i. Contrary to the claim of lB99l the formula of lB92l is not exact 
in the limit of an empty universe. On the contrary, it completely 
fails in this regime: the exact relation in the SCM is then /~^V • 
V = (5, for an arbitrary 5. In fact, this is a general result of dynamics 
in a low-density universe. 

ii. Although the formula of lB92l fails for ^ 0, where it was 
expected to work best, for realistic values of SI (say, SI > 0.1), it 
describes very well the SCM velocity-density relation in voids. It 
also works for overdensities up to 5 ~ 2 - 3. 

The velocity-density relation in the SCM is given in a para- 
metric form. Our goal here was to eliminate this parameter (at least 
approximately) and to provide the relation analytically. We aimed 
at describing the relation in the whole range p G (0, 00) (realisti- 
cally, up to pvir). Therefore, instead of expanding it around p = ph, 
we adopted an entirely different approach. Namely, we derived 
asymptotes of the relation in the highly non-linear regime: p/pb 3> 
1 (5 > 1) for overdensities and pb/p > 1 (0 < 1 + 5 < 1) for 
voids. (Although we also 'expanded' around them, in a sense that 
we also calculated next-to leading-order terms.) These two asymp- 
totes turned out to be qualitatively different. Inspired by their func- 
tional forms, we invented semi-phenomenological fits to the exact 
relation (separately for overdensities and voids), fulfilling the linear 
theory condition /^^V • v = 5. 

For overdensities, our main result is Formula ( 15 lb . It describes 
well the exact relation in the SCM up to the turn-around (for 
O, = 0.25, 5ta = 13.5). As already stated, the formula of IB92| 
starts to deviate from the exact relation for 5 ~ 3. We have also fit- 
ted the regime 5 G {Sta., 5vir), though virialisation and departures 
from spherical symmetry make practical applicability of the SCM 
in this regime very limited. 

In case of voids, the most important results of this paper are 
Formulae ( I41b and J44b . with Gmin given by Equation ( 142b . Com- 
pared with the SCM, simple Formula J41b has a maximal error of 
about 2% and is probably sufficient for practical applications. The 
formula of iB92i is an even better approximation, except for the limit 
(5-^—1, where for Q, — 0.25 it has approximately 5% relative er- 
ror. Our more complicated formula ( 144b is extremely accurate in 
the whole range S < 0: its maximal error is about 0.2%. 

An ultimate goal of studies such as the present one is to find 
the relation valid for reali stic r andom cosmic velocity and density 
fields. Unlike the work of lB92l our calculations were greatly sim- 
plified by the strong assumption of spherical symmetry. There is 
therefore no guarantee that better agreement with the SCM implies 
better agreement with the real relation. In order to check this is- 
sue we compared our formulae to fits to results of cosmological 
numerical simulations, that are present in the literature. We have 
found that in voids, our formulae, as well as the formula of lB92l 
describe well the real relation. This is partly a consequence of the 
fact that voids are more spheric al tha n overdensities. In overdensi- 
ties, both our formula and that of'B92' require modification, but ours 
less. This discrepancy is not a failure of the latter of the two, since it 
has never been intended to work for S^2. Our formula J51b . hav- 
ing the same second-order expansion as the formula of iB92i . can 
be regarded as its extension into the mildly non-linear regime (for 
5 up to the turn-around). We have also discussed how to (slightly) 
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modify our formula to better fit numerical simulations. 

In Section[T]we have enlisted arguments for weak dependence 
of the velocity-density relation on the cosmological parameters. 
Therefore, in the present analysis we set A = 0. To study the limit 
^ we have retained fi-dependence of the equations of the 
SCM. Analysing these equations we have confirmed that for real- 
istic values of Q, the S7-dependence of the relation is indeed very 
weak. In final formulae it has been therefore neglected, except for 
Formula l l42t for Gmin. The difference between 0min for = 1 
and Q = 0.25 is about 5%. In fact, if we want to have better accu- 
racy, there is no guarantee that A does not contribute at a compara- 
ble level. It is then worth to repeat the analysis with A = 1 — il. 
We plan to undertake such a study in the future. 
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APPENDIX A: THE VELOCITY-DENSITY RELATION IN 
AN EMPTY UNIVERSE 

The general equation of motion for the cosmic pressureless fluid in 
comoving coordinates is 



9v 1 , a 

— + - (v ■ V) v + -V = g, 

ot a a 

where g is the peculiar gravitational acceleration, 

, , ^ f ,3 , 5(x',t)(x'-x) 



(Al) 



(A2) 



(e.g. Peebles 1980). For \S\ <C 1 we can neglect the non-linear 
term on the LHS of Equation jAU . Let us choose some instant 
of time, ti, of the evolution of an open universe when already 
f7 <C; 1. For such Q, perturbations stop growing, so for t > ti, 
g(x, t) = gi(x)/a^. Our Equation dAIb simplifies then to 



-[av(x,t)] = ^. 
The solution is 



(A3) 



a(t) a(t) a(t) 

where the conformal time rj is in general defined in Equation 
The last term in Equation l lA4b is the homogeneous part. Here we 
do not assume a priori irrotationality of the velocity field, so we 
retain this term. (Though it does not contribute to the velocity di- 
vergence, because V ■ F = 0.) The limit Q corresponds to 
Tj ^ oo. Therefore, in the above equation we can neglect the terms 
Vi/a and F/a, as well as rji. This yields 

■v^Ho^riag. (A5) 
From Equation l lA2b we have 



v,(x) ^ F(x) 



(A4) 



V • g = -AnGpboS = -^H^naS . 
This yields in iA5i 



(A6) 



(A7) 



In an (almost) empty universe H{t) = and a{t) — t/to, hence 
Ha — Ho. Also, the general relation l|5) between Q and the con- 
formal time simplifies then to = 4e~^ . Substituting this in Equa- 
tion ( IA7t we obtain V-v = — iifo6r;e~''5. Comparing this equation 



Velocity-density relation 1 1 

with Equation ([25}, we identify the factor Qr]e^^ as the low-Sl limit 
of the factor f{0.). Hence, 

V-v = -f/o/(J^)5, (A8) 

in agreement with the general linear theory prediction. Equa- 
tion Q. 

Now, we claim that in the limit SI — » 0, solution iA5l is also 
a solution of the general equation of motion dAlt . for arbitrary S. 
To prove this statement we have to demonstrate that in this limit, 
the non-linear term in equation dAlt is negligible. Substituting so- 
lution JA5b in this term gives 

1^ + -V = g - i (i/,7Sag • V) (Ho" Sag) ■ (A9) 

The amplitude of the second term on the RHS is of order 
H^^rfa y V • g ~ HQ^ffa g H^Qa S. The amplitude of the sec- 
ond term relative to the first is thus 

^ ~ Ho^iHafri^n 5 rfQ.5 ^ T^e-^^S , (AlO) 

and in the limit f2 ^ it tends to zero. (Formally speaking, for 
arbitrary e > and arbitrary S, there always exists 77^ such that for 
all ri > rjc, ri^e~'^\5\ < e.) Thus, in the limit Q —* the non-linear 
term in the equation of motion becomes negligible, for arbitrary 
value of S. This is why in every matter-only, open universe, the 
velocity-density relation evolves towards the linear one. 



APPENDIX B: THIRD-ORDER EXPANSION FOR IN 
VOIDS 

Our aim here is to extend calculations of Section |5] for voids up 
to third order in the perturbation parameter (p (4> is assumed to be 
large, but not infinitely large). We begin applying to Equation dllt 
the equality sinhi^ = cosh<j) — exp(— 0) and expand this equa- 
tion up to terms of the order cosh~^ (f>. Solving perturbatively the 
resulting equation for (j)z we obtain 

v,w, 3.^^ - 1001 + 10 

cosh 03 = cosh <^2 — ] , (Bl) 

cosh 01 

where cosh <j)2 is given by Equation ( I37t . 4>i by d39l l and the second 
term on the RHS of the above equation is a small correction. This 
enables us to write 

Vl + cosh 03 ^ + cosh - " ^^f' + ^° . (B2) 

2 cosh"''"' <j>i 

Using the above equation in Equation dl7t yields 

3 In {sinhri-rif 



{l + 5f, (B3) 

or, finally. 
Here, 

F{5,ri) = 31n^(l +(5) + [10- 61n(2p)] ln(l + 5) + 31n^(2ff) 
- 10 ln(2p) -no , (B5) 

and 51(77) is given by Equation d36t . Inspecting terms in the above 
equation we see that Equation l lB4t can be indeed written in the 
form (gill. 



